RATE OF CONVERGENCE OF THE NANBU PARTICLE SYSTEM FOR 

HARD POTENTIALS 



NICOLAS FOURNIER AND STEPHANE MISCHLER 



Abstract. We consider the (numerically motivated) Nanbu stochastic particle system associ- 
ated to the spatially homogeneous Boltzmann equation for true hard potentials. We establish 
a rate of propagation of chaos of the particle system to the unique solution of the Boltzmann 
equation. More precisely, we estimate the expectation of the squared Wasserstein distance with 
quadratic cost between the empirical measure of the particle system and the solution. The rate 
we obtain is almost optimal as a function of the number of particles but is not uniform in time. 



1. Introduction and main results 

1.1. The Boltzmann equation. The Boltzmann equation predicts that the density f{t,v) of 
particles with velocity u G at time t > in a spatially homogeneous dilute gas solves 

(1.1) dtft{v)^]-[ dvJ doB{\v~vA,0)[ft{v')ft{v:)- ft{v)ft{v,)\, 



§2 



2 

where the pre-coUisional velocities are given by 



(1.2) u = w (w,w^,cr) = — — + CT, v^^v^{v,v^,(j) = — a 

and 9 — 9{v,v„,a) is the deviation angle defined by cos 6 = • a. The collision kernel 

B{\v — v„\,9) > depends on the nature of the interactions between particles. See Cercignani 
[TT] , Desvillettes [H] , Villani jJT] and Alexandre T for physical and mathematical reviews on this 
equation. Conservation of mass, momentum and kinetic energy hold at least formally for solutions 
to (jl.ip and we classically may assume without loss of generality that J^^^ fQ{v)dv — 1. 

We will assume that the collision kernel is of the form 

(1.3) B{\v~-v4,9) sin9 = <P{\v-v^\)Pi0) with /3 > on (0,7r/2) and /3 = on [n/2,7:]. 

This last condition (3 — on (7r/2,7r] is not a restriction, since one can always reduce to this case 
for symmetry reasons, as noted in the introduction of Alexandre et al. [3'. 

When particles behave like hard spheres, it holds that ^{z) — z and /3 = 1. When particles 
interact through a repulsive force in l/r*, with s e (2, cxd), one has 

$(z) = z^i with 7 = ^— ^ e (-3, 1) and m) - zsiQ-^-" with v = e (0,2). 

s — 1 s — 1 

One classically names hard potentials the case when 7 e (0,1) (i.e., s > 5 and ly e (0,1/2)), 
Maxwell molecules the case when 7 = (i.e., s — 5 and i> = 1/2) and soft potentials the case when 
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7 G (—3, 0) {i.e., s G (2, 5) and v € (1/2, 2)). The present paper concerns Maxwell molecules, hard 
potentials as well as hard spheres, so that we always assume 7 G [0, 1]. 

1.2. Stochastic particle systems. As a step to the rigorous derivation of the Boltzmann equa- 
tion, Kac [26] proposed to show the convergence of a stochastic particle system to the solution 
to (jl.ip . Kac's particle system is a (R'^)^-valued Markov process with infinitesimal generator Cn 
defined, for </> : (M^)^ R sufhciently regular and v = (wi, . . . , wat) G (M^)^, by 



For h G R'^, we note hei = (0, . . . , 0, /i, 0, . . . , 0) G (R'^)^ with h at the i-th place. Roughly speaking, 
the system is constituted of N particles entirely characterized by their velocities {vi, . . . ,vn) and 
each couple of particles with velocities {vi,Vj) are modified, for each a G at rate B{\vi — 
Vj\,0)/ {2{N — 1)) and are then replaced by particles with velocities v' {vi,Vj,a) and v'^{vi,Vj,a). 

In the present paper, we will consider a slightly modified and non-symmetric particle system 
introduced by Nanbu [34] . The Nanbu stochastic particle system corresponds to the generator Cn 
defined, for cj) : (M^)^ ^ R sufhciently regular and v = (wi, . . . , vn) G (M^)^, by 




'i - Vj\,6)da. 



This system still describes N particles characterized by their velocities {vi, . . . ,vn), but now each 
couple of particles with velocities {vi,Vj) are modified, for each a G at rate B{\vi — Vj\,9)/N 
and are then replaced by particles with velocities v'{vi,Vj,a) and Vj. Thus only one particle is 
modified at each "collision", but the rate of collision is multiplied by 2. All in all, the asymptotic 
behavior, a.s N ^ 00, should be the same. 

1.3. Aims. Our aim is to prove that as N tends to 00, the Nanbu stochastic system is asymptoti- 
cally constituted of independent particles with identical law governed by the Boltzmann equation, 
and better, to quantify this convergence. 

There are two main motivations for such a study, (i) From a physical point of view, we want to 
know how well the Boltzmann equation approximates true particles. Of course, true particles are 
subjected to classical (non random) dynamics. However, Kac's particle system can be seen as such 
a (strongly) simplified system. In some sense, the Nanbu non-symmetric particle system under 
study is again a simplification of Kac's dynamics, (ii) From a numerical point of view, we want to 
know how well the particle system approximates the Boltzmann equation. It is then important to 
get rates of convergence, to know how to choose the number of particles (and the cutoff parameter) 
to reach a given accuracy. 

One might be surprised at first glance: why approximate a Boltzmann equation by a particle 
system which the Boltzmann equation is expected to approximate? The reason is of numerical 
computation order. Indeed, it is numerically interesting to approximate a particle system with a 
huge number of particles (as in a true gas) by a particle system with much lesser particles. 

The main difficulty lies in the fact that even if the particle system is initially constituted of 
independent particles, they do not remain independent for later times, because of interactions. 
Hence to answer the convergence issue, we have to prove that particles asymptotically become 
independent and in the same time to identify their common law: we have to prove that the system 
is chaotic in the sense of Kac [55] . 
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We are able to prove and quantify the chaotic property for Nanbu's particle system. Unfor- 
tunately, our study does really not seem to work for Kac's particle system. From the physical 
point of view, Nanbu's system is less pertinent. However, we believe that the behaviors of the two 
systems are very similar, so that our results should also hold true for Kac's particle system. From 
the numerical point of view, both systems are expected to approximate the solution to (jl.ip with 
an error of the same order, so that the system under study is as interesting as Kac's system. 

We will also study a cutoff version of Nanbu's system. This is motivated by two reasons. From 
a numerical point of view, the particle system with generator Cm cannot be directly simulated, 
because each particle collides with infinitely many others on each time interval (except for hard 
spheres). Thus we have to introduce a cutoff. From a technical point of view, we are not able 
to prove directly our estimates for the particle system without cutoff: we have to study first the 
particle system with cutoff and then to pass to the limit. 



1.4. Assumptions. We assume that the collision kernel is of the form (|1.3|) with 

(1.5) 3 7 e [0, 1], Vz > 0, $(z) = z^, 
and either 

(1.6) V0e(O,V2), /3(^?) = 1 
or 

(1.7) 3 e (0, 1), 3 < Co < ci, Ve* e (0,7r/2), cqQ-^-" < I3{6) < cid-^-" . 

This work could probably be extended to G (0, 2), since the important computations on which 
it relies also hold in this case. However, this would introduce several technical difficulties. Since 
Maxwell molecules and hard potentials, which we study, satisfy (|1.7p with v G (0, 1), we decided 
to avoid these technical complications. 

The propagation of exponential moments requires the following additional condition 

(1.8) j5{9) — b{cosd) with b non-decreasing, convex and on [0, 1). 

In practice, all these assumptions are satisfied for Maxwell molecules (7 = and ly ~ 1/2), hard 
potentials (7 e (0, 1) and ly e (1, 1/2)) and hard spheres (7 = 1 and /3 = 1). 

1.5. Notation. For 9 G (0,7r/2) and z e [0, 00) we introduce 

(1.9) H{9)= / I3{x)dx and G{z) = H-^{z). 

Je 

Under (|1.7p . H is a continuous decreasing bijection from (0,7r/2) into (0, 00), and its inverse 
function G : (0,oo) 1-^ (0,7r/2] is defined by G{H{d)) = 0, and H{G{z)) = z. It is immediately 
checked that under p.7p . there are some constants < C2 < C3 such that 

(1.10) Vz>0, C2(l + z)-i/'' <G(z) <C3(l + z)-i/'^ 

and, as checked in [191 Lemma 1.1], there is a constant C4 > such that for all x, y G M+, 

(1.11) r {G{z/x) - G{z/y)f dz < C4 ~ . 

J a x + y 

Under (|1.6|) . we have G{z) = (7r/2 — z)_|- (with the common notation x^ — maxjx, 0}) and a direct 
computation shows that (|l.lip also holds true. 
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1.6. Well-posedness. Let VkO^^) be the set of all probability measures / on M.^ such that 
/r3 {vl'^fidv) < oo. We first recall known well-posedness results for the Boltzmann equation, 
as well as some properties of solutions we will need. A precise definition of weak solutions is stated 
in the next section. 

Theorem 1.1. Assume (dH]) and or dTT]). Let /o G p2(R^)- 

(i) If J ~ 0, there exists a unique weak solution {ft)t>o G C([0, oo), 'P2(M'^)) to (jl.ip . // 
/o e Vp{R^) for some p > 2, then supfg^oo) /r3 \v\^Mdv) < oo. If J^^ fo{v) log fo{v)dv < oo 
o'' */ /o G ^^{MP) and is not a Dirac mass, then ft has a density for all t > 0. 

(ii) If "/ £ (0, 1] , assume additionally (jl.Sp and that 

(1.12) 3pe(7,2), / el''l7o(d^^) < oo. 

There is a unique weak solution {ft)t>o G ^^([0, oo), 'P2(IR'^)) to p.ip such that 

(1.13) Vqe(0,p), sup [ el"IVt(dw) < oo. 

[0,oo) JR3 

Under (jl.7p anrf i/ /o is rtof a Dirac mass, then ft has a density for all t > 0. Under (jl.6p and if 
/o /las a density, then ft has a density for all t > 0. 

Concerning well-posedness, see Toscani-Villani ^] for Maxwell molecules, [22j [M] for hard 
potentials and 5, 33, 271 lEl [2H] for hard spheres. The propagation of moments in the Maxwell 
case in standard, see e.g. Villani [51] Theorem 1 p 74]. The propagation of exponential moments 
for hard potentials and hard spheres, initiated by Bobylev [J, is checked in [22l[28]. Finally, the 
existence of a density for ft has been proved in 17 (under (|1.7p and when /o is not a Dirac 
mass and belongs to 7^4 (R'^)), in [33 (under (II. 6p when /o has a density) and is very classical by 
monotonicity of the entropy when /o has a finite entropy, see e.g. Arkeryd [3]. 

We now introduce our particle system with cutoff. 

Proposition 1.2. Assume ([T3|) and pTE)) or pTT)) . Lei /o G 7'2(K^) a^rf a number of 

particles N > 1 be fixed. Let (Vjj')i=i,...Ar be i.i.d. with common law /q. 

(i) For each cutoff parameter K G [l,oo), there exists a unique (in law) Markov process 
{Vt'^'^)i=i^,,,^N,t>Q with values in (ffi^)^, starting from (Vjj*)i=i^...jv and with generator Cn,k de- 
fined, for all bounded measurable </> : (R'^)^ H> R and any v = (ui, . . . , vn) G R'^, by 

CN,K(f>{v) = T^y^ / [4>{^ + {v'{vi,Vj,a) -Vi)ei) - (j>{v)]B{\vi - Vj\,e)l^g>c{K/\vi^vi\-')}dcr, 

with G defined by (|1.9p anrf, /or /i G R^, /le^ = (0, . . . , /i, . . . , 0) G (R^)^ with h at the i-th place. 

(ii) There exists a unique (in law) Markov process (V^''^'°°)i=i,...,Ar,t>o with values in (M.^)^ , 
starting from (V^*)i=i,...jv Cind with generator Cn defined, for all Lipschitz bounded function <f) : 
(R3)^ H>R and any v = {vi,...,vn) G R^, by ((TTi)) . 

Let us emphasize that the cut-off used for defining the generator Cn,k is not the usual one since 
it depends not only on the deviation angle 9 G (0, 2it) but also of the relative velocity |u — It is 
more convenient in order to perform the computations we want to do. It might also be convenient 
for practical simulations. Indeed, the total rate of collision of the particle system does not depend 
on the configuration of the velocities: it always equals 2tt{N — 1)K. Hence, the (mean) simulation 
cost of the particle system on a time interval [0,r] is proportional to {N — 1)KT. 
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1.7. Wasserstein distance. For g,g ^ V2{^^), let ^{{gjg) be the set of probability measures on 
M.^ X M.^ with first marginal g and second marginal g. We then set 

W2{g,~g) = infjf/ \v-vfT^idv,dv))^ ; 6 ?^(.9,5)| • 

[ \JR3xR3 / J 

This is the Wasserstein distance with quadratic cost. It is well-known that the inf is reached. 
We refer to Villani Chapter 2] for more details on this distance. A remarkable result, due to 
Tanaka [351131]: is that in the case of Maxwell molecules, 1 1— W2(/t, ft) is non-increasing for each 
pair of reasonable solutions /, / to the Boltzmann equation. The present work is strongly inspired 
by the ideas of Tanaka. 

1.8. Empirical law of large numbers. For / £ ■P2(K'^) and > 1, we define 



(1.14) eAr(/):=E 



Wl I f,N 



N \1 



with Xi, . . . , Xn independent and /-distributed. 



Since a /-chaotic stochastic particle system is asymptotically constituted of i.i.d. /-distributed 
particles, £jv(/) is the best rate (as far as >V| is concerned) we can hope for such a system. We 
summarize in the following statement the best available estimates on ejv(/), essentially picked up 
from Rachev-Ruschendorf [351 Theorem 10.2.1], [31] Lemma 4.2] and Boissard-Le Gouic [51. 



Theorem 1.3. For all A > 0, all k>2, all f e Pfe(M^) verifying J^^ \v\''f{dv) < A, all N >1, 

(1.15) ENif) < 



CA,kN-'-' forr = {k~2)/{bk) z/fce(2,7), 

CA,k,rN-^'' for any r < 1/7 if k ^ 7, 

C^V^-^r for r = 1/7 if ke {7,20], 

CA,k,rN-^''' for any r < {k -2)/{6k + 6) if k > 20. 



The proof of Theorem II. 31 will be discussed in Appendix \K\ We shall use different proofs for the 
different cases, which explains why the statement is awful. The above rates are clearly not optimal. 
If / has infinitely many moments, then we almost obtain (/) < CN-'^^^. This power 1/3 might 
be optimal. See Barthe-Bordenave \W and the references therein for a discussion concerning this 
issue. We also refer to [25, Theorem 2.13] (and the remarks which follow) for a general discussion 
about the rate of chaoticity for independent and dependent random arrays. 

1.9. Main result. Our study concerns both the particle systems with and without cutoff. It 
is worth to notice that for true Maxwell molecules and hard potentials, G (0,1/2] so that 
I — 2/1/ < —3 and the contribution of the cut-off approximation vanishes rapidly in the limit 
AT 00. 

Theorem 1.4. Let B be a collision kernel satisfying (|1.3p . p.Sp and (|1.6p or p.7p and let fo £ 

^2(1^"^) riot be a Dirac mass. If ^ > 0, assume additionally (|1.8p and (jl.l2l) . Consider the unique 
weak solution {ft)t>o to (|l.ip defined in Theorem \1.1\ and, for each N > 1, K Cz [l,oo], the unique 
Markov process {V/'^'^)i=i^,,,^N^t>o defined in Proposition ll.SX Let fi^'^ := N^^ 5yi,N,K . 

(i) Maxwell molecules. Assume that 7 0, (|1.7|) and either J^^ fo{v)log fo{v)dv < 00 or 
fa G P4(R'^). There is a constant C such that for all T > 0, all N > 1, all K G [1, 00], 

(1.16) sup E[>V|(/if ft)] <C(1 + r)2 sup SNift) + CTK^-^^r 

[0,T] [Q,T] 
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If fo e Pk{R^) for some k > 2, we have supjQ j^-j J^^ \v\''ft{dv) < oo and we can use Theorem ] l.S\ 
to bound sup[g £N{ft)- In. particular if fo £ Vki^'^) for all k > 2, then for all e > 0, all T > 0, 
all N >1, all K e [l,oo], 

(1.17) sup E[Wi{fi^'^, ft)] <Ce(l + T)27V-i/3+e ^ CTK^-'^/". 

[0,T] 

(ii) Hard potentials. Assume that 7 G (0, 1) and (|1.7p . For all e G (0, 1), all T > 0, there is a 
constant Cs^^t such that for all N > \, all K 00], 

(1.18) sup E[W2'(//f ft)] <C,,T f sup ENift) + K^-^''^ . 

[0,T] \[Q.T\ j 

Consequently, for all e G (0, 1), all T > 0, there is C^^t such that for all N >1, all K ^ [1, 00], 

(1.19) sup E[>V|(/i^'', /*)] <Ce,T(iV-l/3 + Xl-2/-)l-. 

[o,r] 

(Hi) Hard spheres. Assume finally that 7 = 1, (jl.6p and that fo has a density. For all e G (0, 1), 
all T > 0, all q (z (l,p), there is a constant C^^q^T such that for all N > 1, all K d [1, 00), 

(1.20) sup E[Wi{fi^''',ft)]<C,.q,T((snpeN{ ft)] + e"^' I e^--^^ 

[O.T] I V [0,T] / 



Thus for all e G (0,1), all T > Q, all q ^ (1?^)? there is Ce,q,T such that for all N > 1, all 

Ke[i,cx^), 

(1.21) supEiWK/if''',/*)] <a,,,T(iV-^/=^+^ + e-^')e^---^. 

[0,T] 

Concerning the rate of convergence of the simulation algorithm, we have the following. 

Remark 1.5. Recall that the simulation cost per unit of time is proportional to {N — 1)K . 

(i) For Maxwell molecules and hard potentials the error (forW2) is{N-^/^+K^/^-'^/'')^- . For 
a given simulation cost t, the best choices are N ~ t-(6-3i^)/(6-2!/) ^ .^iy/(6-3u) ^ which leads 
to an error in t--(2-i^)/(12-4!^)+^ true hard potentials and Maxwell molecules, this is at worst 
T-3/20+ and at bestT-^'''+. 

(ii) For hard spheres, make the choice K ~ (logTV)" with a G (1/(7,1). Then e^^ << N'^ 
for any e G (0, 1) and e"^ « iV^'' for any r > 1. With this choice, we thus find an error in 
^-i/6+£ j-pj. ^ simulation cost in N{\ogN)'^. Consequently, for a given simulation cost t, we find 
an error in t^^/^+. 

We excluded the case where fo is a Dirac mass because we need that ft has a density and 
because if fo = S^^ , then the unique solution to (11.11) is given by ft = S^g and the Markov process 
of Proposition 11.21 is nothing but v/'^''^ — {vo,...vo) (for any value of K G [l,oo]), so that 

' = dyg and thus W2ift,fJ't ' ) = 0. 

1.10. Comments. We thus show that the empirical law of the particle system converges to ft as 
fast as i.i.d. /(-distributed particles (up to an arbitrary small loss if 7 7^ 0). This is thus almost 
optimal in some sense. However, this is optimal only as far as W2 is concerned: we would have 
preferred to work with another distance and to obtain a rate in iV~^/^ as is expected for laws of 
large numbers. Here we obtain a rate in N~^/^ , since W2 is squared. However, W2 enjoys several 
properties that make it quite convenient when studying the Boltzmann equation, mainly because 
of the role of the kinetic energy. Another default of this work is that we obtain a non-uniform 
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(in time) bound. For Maxwell molecules, the bound is slowly increasing (as T^) but for hard 
potentials, it is growing very fast. 

Note also that for hard spheres, we are not able to treat the case where K = oo: we need to 
let K and N go to infinity simultaneously, with some constraints. We believe that this is only a 
technical problem, but we were not able to solve it. However, we still obtain a very reasonable rate 
of convergence (as a function of the computational cost). 

Our proof is based on a coupling argument: we couple the iV-particle system with a family 
of N i.i.d. Boltzmann processes, in such a way that they remain as close as possible. We prove 
an accurate control on the increment of the distance between the two systems at each collision. 
This last computation is similar to those of [1H[T^ concerning uniqueness of the solution to 
However, we need to handle much more precise computations: in |22| . when studying the distance 
between two solutions to (II. both were supposed to have exponential moments. Such exponential 
moments are known to propagate for solutions to since the seminal work of Bobylev [7] , but 
for the particle system under study, we are not even able to prove the finiteness of a moment of 
order 2 + £, e > 0! We thus need a very precise refinement of the computations of [2^ [T^ . 

All these problems do not appear when studying Maxwell molecules. Roughly, the collision op- 
erator is globally Lipschitz continuous for Maxwell molecules and only locally Lipschitz continuous 
for hard potentials (which explains why large velocities have to be controlled by using exponential 
moments). This is why we obtain a better result for Maxwell molecules. 

Note that for the (physically more relevant) Kac particle system moments are known to prop- 
agate (uniformly in N), see Sznitman [3B] and also [3T], which would simplify greatly the proof 
at many places. However, we are not able to exhibit a suitable coupling. This is due to the fact 
that in Kac's system, each collision modifies the velocity of two particles. In Nanbu's system, the 
Poisson measures governing two different particles are independent, which is not the case for Kac's 
system (because each time a particle's velocity is modified, another one has to be also modified) 
although the larger is the number of particles, the lower the correlation is. As a consequence, it is 
more difficult to couple the A^-particle symmetric Kac's system with N independent copies of the 
Boltzmann process and we did not succeed. 

1.11. Known results. Such a chaos result for the Boltzmann equation with bounded cross section, 
or for related models, has been first established without any rate by Kac [25 (for the so-called 
Maxwell molecules Kac's model which is roughly a "toy one-dimensional" Boltzmann equation) 
and then by McKean [S^ and Griinbaum 24 . For unbounded cross section, the chaos property 
has been proved by Sznitman |36j for hard spheres, still without rate. 

For Maxwell molecules with Grad's cutoff, a nice rate of convergence (of order 1/N in total 
variation distance on the two-marginal) has been obtained by McKean [2 9) and improved by 
Graham-Meleard [33]. This was extended by Desvillettes-Graham-Meleard [13], see also ^U\, to 
true (without Grad's cutoff) Maxwell molecules, but with a rate in N~^e^^ + (-^^ith 
the notation of the present paper). From a numerical point of view, this leads to a logarithmic 
convergence as a function of the computational cost. 

More recently, a uniform in time rate of chaos convergence of Kac's stochastic particle system to 
the Boltzmann equation for two unbounded models has been established in [31] [10] (see also [32]). 
by taking up again and improving Griinbaum's approach. For true Maxwell molecules, uniform in 
time rate of convergence of order N''^/ ^ for any (5 > 0, for a weak distance on the two-marginals 
has been proved in 13 H Theorem 5.1] when the initial condition /o has a compact support. This 
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result was improved and made more precise in |10t Step 3 of the proof of Theorem 8], where, still 
for true Maxwell molecules, uniform in time rate of convergence of order iV"^/^''^, for the same W2 
Wasserstein distance as used in (|1.16|) . has been proved for any initial condition /o satisfying (|1.12p . 
Hard spheres have also been studied in [311 Theorem 6.1]; a uniform in time rate of convergence 
of order l/(logA^)" with a > small, for the Wi distance on the two-marginals has been proved. 
When applying the methods of [31] [351 [TUj on finite time intervals, the previous rates can not be 
really improved. Finally, let us mention that the present work follows some of the ideas of [18] . 
which concerns the Kac equation. 

To summarize: 

• We obtain the first rate of convergence for hard potentials and this rate is reasonable. Recall 
that hard potentials are twice unbounded (the velocity cross section is unbounded and the angular 
cross section is non-integrable) , while Maxwell molecules enjoy a bounded velocity cross section 
and hard spheres an integrable angular cross section. 

• For hard spheres and Maxwell molecules, we prove a much faster convergence than [3T1I3^[TU] . 
but we are restricted to finite time-intervals and we cannot study Kac's system. 

Let us finally mention that we use a coupling method, as is widely used since the famous cours 
a I'ecole d'ete de Saint- Flour by Sznitman |37) for providing rate of chaos convergence for the so- 
called McKean-Vlasov model and that such methods have been recently adapted to non-globally 
Lipschitz coefficients by BoUey-Cahizo-Carrillo in [5], making use of exponential moments. 

1.12. Plan of the paper. In Section[2l we make precise the notion of weak solutions, rewrite the 
collision operators in a suitable form and check a accurate version of a lemma due to Tanaka [39] . 
Section [3] is devoted to the cornerstone estimate on the collision integral. In Section |4] we prove 
the convergence of the particle system with cutoff. The cutoff is removed in Section [S] Finally, 
elements of the proof of Theorem 11.31 on quantitative law of large number for empirical measures 
are presented in Appendix [Xj 



2. Preliminaries 



2.1. Rewriting equations. We follow here [30]. For each X G K'^, we introduce I{X), J{X) e R'^ 
such that (j^, l^T' "^W ^ ^ direct orthonormal basis of K'' and, of course, in such a way that 
/, J are measurable functions. For X, w, G R"^, for 9 G (0, 7r/2) and G [0, 27r), we set 

f r(X, (cos <^)/(X) + (sin V3)J(X), 



a{v,v^,9,ip) 



1 — cosO , , sin0 
[v - v^) H —T[v - v^,(p), 



(2.1) 



which is a suitable parametrization of ()1.2|) : write u G S as ct 



^ v'{v,v^,9,ip) V + a{v,v^,6,ip), 



\v—v^ I 



\v—v^ I 



sin 9 cos ip - 



sin sin Let us define, classically, weak solutions to ()1.1|) . 

Definition 2.1. Assume and (HI]) or ([LT]). A family {ft)t>o e C([0, 00), 7^2(1^^)) 

is called a weak solution to (jl.ll) if it preserves momentum and energy, i.e. 



(2.2) 



Vt > 0, 



vft{dv) 



vfo{dv) and 



\v\'.ft{dv) 



\v\^fo{dv) 
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and if for any : R'^ M> M bounded and Lipschitz-continuous, any t G [0,T], 

(2.3) / ^{v)ft{dv)= f ^{v)fo{dv)+ f f f A(f>{v,v.)fs{dv,)fs{dv)ds 



where 

(2.4) Acl){v,v,) = \v-v,p P{d)de dip[(t){v + a{v,v^,e,^p))-4){v)]. 

Jo Jo 

Noting that \a{v,v^,,9,(p)\ < C6\v — w*| and that J^^^ 9l3{9)d6, we easily get \A(j){v,v^)\ < 
C^\v — v^\'^^^ < C^{1 + \v — w*P), so that everything makes sense in (|2.3p . 

We next rewrite the colhsion operator in a way that makes disappear the velocity-dependence 
|u — P in the rate. Such a trick was already used in [3T] and P^ . 

Lemma 2.2. Assume ([Ol) . (IT3)) and or pr7| . Recalling and (pH . define, for 

z e (0, oo), If e [0, 27r), V, G and K e [1, cx)), 

(2.5) c{v,v^., z,(p) := a[v,v^,,G{z/\v — v^,\'^),(p] and ck{v,v^, z, (p) := c{v,v^,, z,ip)TL^z<^j(y 
For any bounded Lipschitz : R'^ i-> M, any v, G 

(2.6) ^0(w,w,) = J dz j dip{4'[v + c{v,v^:, z,ip)] ~ (f)[v\^. 

For any > 1, G [1, oo), v — (wi, . . . ,vn) G (M^)^, any bounded measurable : (R'^)^ M> R, 

(2.7) ^N,K(t)i'y) ^ dz dip[(j){-v + CK{v^,Vj,z,ip)e,) - 

r^j Jo Jo 

For any N > 1, any v ~ [vi, . . . , uat) G (R'^)"'^, any bounded Lipschitz (f) : (JS?)^ ^ M, 

(2.8) ^ ITt^ I / rf</'[0(v + c(wj,Uj,z,93)ej) -(/)(v)]. 

^ Jo Jo 

Proof. To get (12. 6p . start from (|2.4p and use the substitution = G{z/\v — w*!"*") or equivalently 
H{d) = z/jw — which implies Iv — v^^p (3{9)d9 = dz. The expressions (|2.7[) and (|2.8p are checked 
similarly. □ 

2.2. Accurate version of Tanaka's trick. As was already noted by Tanaka [39., it is not possible 
to choose / in such a way that X i— >■ I{X) is continuous. However, he found a way to overcome 
this difficulty, see also [IHl Lemma 2.6]. Here we need the following accurate version of Tanaka's 
trick. 

Lemma 2.3. Recall (j2.ip . There are some measurable functions LpQ,ifi : R'^ x R'' i— )■ [0, 27r), 
such that for all A, F G R^, all ip G [0, 27r), 

r(A, if,) ■ r(y, ^ + ipo{X,Y)) = X ■Ycos'if + fi{X,Y)) + \X\\Y\ sm\f + (^i(A, F)), 

|r(A,^)-r(r,^ + ^o(A,r))| < \x-y\. 

Proof. First observe that the second claim follows from the first one: writing fi — fi{X, Y) 

|r(A, f) - r(r, ^ + ^o)l' =|r(A, if)\^ + |r(y, ^ + ^o)l' - 2r(A, ^) • r(r, ^ + p^) 

= |A|2 + |y|2 - 2(A -Ycos^iif + ipi) + |A||y| siu\f + ^i)) 

<|aP + |yp 2A • r = lA - rp. 
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We next check the first claim. Let thus X and Y be fixed. Observe that T{X, ip) goes (at constant 
speed) all over the circle Cx with radius \X\ lying in the plane orthogonal to X. Let ix € Cx and 
iy G Cy such that X, Y, ix, W belong to the same plane and ix ■ iy = X - Y (there are exactly two 
possible choices for the couple {ixjiy) if X and Y are not collinear, infinitely many otherwise). 
Consider ipx and (py such that ix ■— r(X, ipx) and iy := T{Y, ipy). Define jx '■— r(X, ipx +7r/2) 
and jy := r{Y, (py +tt/2). Then jx and jx are collinear (because both are orthogonal to the plane 
containing X,Y,ix,iy), satisfy jx - jy — IjxWjvl — l-'^H^I and ix - jy — iy-jx = 0. Next, observe 
that T{X, p + px) — ix cos p + jx sin p while T{Y, p + py) = iy cos p + jy sin p. Consequently, 
T{X,p + px) ■ r(r, ^ + py) ^ ix ■ iy cos^ p+ jx ■ jy sin^ p ^ X ■ F cos^ p + \X\\Y\ sin^ p. The 
conclusion follows: choose pq := py — px and pi := —px (all this modulo 27r). □ 



3. Main computations of the paper 
The following estimate is our central argument. 

Lemma 3.1. Recall that G was defined in (II. 9p and that the deviation Junctions c and ck were 
defined in (|2.5p . For any v, u*, w, € M.'^ , any K G [1, oo), 

/ \\v + c{v,v^,, z,p) — V — ck{v,v^,z,p + pq{v — v-^,,v — v.t))\ ~ \v ~ v\'^\dpdz 
<Af (w, w*, w, w,) + A2 {v,v^,v,v^) + (-y,-y,,-D,-D,), 
where, setting <^k(x) = tt /q^^(1 — cos G{z/x"'))dz and ^>Kix) — tt /^(l ~ cos G{z/x'^))dz, 

(?;,«*, =2|u — — w,| / [G'(z/|w — P) — G(z/|{j — Uh-D] dz, 

Jo 

A2 {v,v^,v,v^) = - [{v -v) + (w, - w,)] • [{v - v^)^k{\v - w*|) - (w - v^)^k{\v - I 
{v,v,:,v,v,:) ={\v — + 2|u — v\\v — w*|)5' /<-(!« — w*|). 

Proof. We need to shorten notation. We write x = \v — v^\, x = \v — v^\, pQ — pq{v — v^,v — v-^,), 
c = c{v, v^,z, p), c = c{v, v^, z,p + pq) and ck = ck{v, v*,z, p + po) = cl{^<x}- We start with 



Ak := 



00 /'27T 



"'0 
K r2Tr 



JO 



(^\v + c — V — ck\'^ — \v — v\'^^ dpdz 
(^|cp + _ 2c • £ + 2(u - w) • (c - £)^ dpdz 
(^|cp + 2{v - v) ■ c^dpdz. 



IK Jo 

First, it holds that \c\^ = | - (1 - cos G(z/a;^))(w - v*) + (sin G(z/x''))r(u - v^,p)\'^/4: = 
(1 — cosG(z/a:^))|w — w,|^/2. We used that by definition, see (|2.ip . T{v — v^,p) has the same norm 
as w — and is orthogonal to w — and that (1 — cos 9)^ + (sin^)^ = 2 — 2cos0. Consequently, 
we have 

K r27T pK 

\c\'^dpdz — 'k\v ~ v^.]"^ i {1 — cos G{z/x^))dz — x^^xix). 
Jo 

Similarly, we also have J^^ \c\'^dpdz — a;^$if (i) and J^^ \c\'^dpdz = a;^5'/<-(a;). 
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Next, using that c = —(1 — cos G(z/x'*'))(w — u*)/2 + {smG{z/x^))r{v — v^,,(p)/2 and that 



Jg'^ r{v - v^.,ip)dtp = 0, 



cdipdz — —{v — w*)7r / (1 — cos G{z / x"'))dz — ^{v ~ v^,)^k{x)- 
Jo 

By the same way, J^'^ J^^ cdipdz — —(v — w*)<i>if (i) and cdipdz = — (w — Vt:)"^Kix)- 

Finahy, c • c = [(1 - cos G(z/a;^))(w - v^) - {smG{z/xf j)r{v - v^,^)\ ■ [(1 - cos G(z/£'^))(w - 
v^,) — {smG{z/x''))r{v - v^,ip + ipo)]/4:. Since J^^ r{v - v^,ip)d(p — J^^ r{v — v^,, ip + ipo)dip — 0, 
we get 

/ c • cdip =—(1 — cos G{z/x^)){l — cos G{z/x^)){v — v^) ■ {v — v^) 
Jo 2 

1 /-^^ 

+ -(sinG(z/a;'''))(sinG(z/£^)) / T{v — v^,,ip) ■ T{v - v^,,ip + (po)d'f. 
4 Jo 

Recalling Lemma 12.31 and using that J^^ cos^{(p + fi)dip ~ J^^ sin^(93 + (pi)d(f = tt, we obtain 

r27r _ 
TT 

c • c(i(^ "2" ^ ^ COS G(z/2;''))(w — v^) ■ {v — v^) 

+ ^(sinG(z/x'''))(sinG(z/i''')) [{v — v^) ■ {v — w,) + |f — - -D*]] . 
But G takes values in (0,7r/2), so that, since \v — Vf,\\v — {;,| > (w — w*) • {v — -D,), 

c • cdcp 

>|[(1 _ cosG(z/a;''))(l - cosG(z/i^)) + (sinG(z/a:'^))(sinG(z/iT))](w - w*) • {i ~ i^) 
= - cosG(z/a;'^)) + (1 cos G(z/i''))](f - v,) ■ {i - v,) 
_ |(1 - cos(G(z/x'') - G{z/i^))){v - V,) ■ - V,). 



Using that 7r(l — cos^) < 29 , we thus get 

c-cd'fdz>{v-v^)-{v-v^)^^^^^^-^^^-xx ('^ {G{z/x-')-G{z/i^)fdz. 
Jo Jo 2 Jo 

All in all, we find 

<x^'^k{x) + x^'^k{x) - (u - w,) • (u - ■5,)[$K(a;) + 
+ 2(i; - w) • [(5 - '5,)'I>K(i) - (v - u,)'I>a:(x)] 

+ 2xi / (G(z/a;^) - G(z/i'^))2dz 
Jo 

+ x^^k{x) - 2{v -v) ■ {v - v,)^'if(a:). 

Recalling that x = |w — x = |u — u*!, we realize that the third line is nothing but Af'{v, u*, v, u*) 
while the fourth one is bounded from above by {v,Vt,v,v^:). To conclude, it suffices to note 
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that the sum of the terms on the two first fines equals 

— {v — u*) • [{v — v^,) — {v — -C,) — 2{v — v)]^k{x) 

+ {v - -5,) • [{v — w,) — {v - w,) + 2{v - v)]^k{x) 
= - {v - v^) ■ {{v -v) + (u* - v*))^Kix) + {v - w,) • {{v ~v) + (w, - 



Next, we study each term found in the previous inequafity. We start with the MaxweU case. 

Lemma 3.2. Assume (|1.3p . (jl.Sp with 7 = 0, (jl.7p and adopt the notation of Lemma \3.1\ For 

all K G [1, 00), all v, w*, w, w* 6 K^, 
(i) Af (u, u*, u, u,) = 0, 

(a) A2 {v,v^:,v,v^,) — (k[—\v — + |w* — where (k ~ Jq (1 — cosG(z))dz, 
(^zm; A^iv,v,,v,v,) < C(|t;|2 + |i;,|2 + |5|2)i^i-2/>'. 

Proof. Point (i) is obvious. Point (ii) immediately follows from the fact that 'i'kix) ~ Ck does not 
depend on x. Point (iii) holds true because ^i<-(a;) = tt/^ (1 — cosG{z))dz < tt G'^{z)dz < 
CK^-^'" by (fTTOl) . □ 

The case of hard potentials is much more complicated. The following result gives a possible and 
useful upper bound on the Af- functions. 

Lemma 3.3. Assume (|1.3p . (|1.5p with 7 G (0, 1), (|1.7p and adopt the notation of Lemma \3.1\ 

(i) For all q > 0, there is Cq > such that for all M > 1, all K G [l,oo), all v,v„,v,v„ G K'^, 



This lemma is very technical. The reason is the following. The solution {ft)t>o has bounded 
exponential moments while, on the contrary, the particle system has only a bounded energy (mo- 
ment of order 2). li K d [l,oo), the particle system has all moments finite, which makes all the 
computations licit, but the moments of order strictly greater than 2 are not uniformly bounded 
with respect to K (at least, we were not able to show it). We will use the previous estimates with 
w,^, (and z,) taken from the solution ft and v,v,, taken in the particle system. Thus, it is very 
important that these estimates do not involve powers greater than 2 of v, v^. For example in point 
(i), only appear in the exponential and this is crucial. 

Proof. Using (|l.lip and that Ix"^ — y''\ < 2\x — y\/{x^^'^ + y^^^), we get 



which is A2 



{v, w,, V, w,) as desired. 



□ 




(3.1) 




{v, w*, w, V*) <2c4\v — — 5,1 



<8C4 



{\v-v^\V\v-v4y—' 



u — w, 1 A I'D — ■0*1 
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Now for any Af > 1, this is bounded from above by 

^i\v -v*\-\v - v^f + Sc4{\v - w*| V lu - V:,\)'^^^TL,^ |.-.,|A|€,-ii.i ^a/-, 

2 <-*^*(|„-„.|v|S-i.|)i-T^ 2 i 



M ,. , ~ u2 o 

< — (|w - w| + |w* — + 8C4 



{\v — w,| A |u — 



2 + 7 



1-T 

If k-D.lAlii-i),! ^ Af . 
^ (|^-^,|V|i-i,|)i-T - 16c4 J 



2 + 7 



M 

<Af(|t; - + Iw, - + 8c4 [16c4(|w - A |w - a{(|,_„.|/,|5_^.|)-,>_M_j 
<Af(|v - + - + 8c4 [16c4(|v| + a{(|„|+|„^|)7>^} 
Fix now q > and observe that 

Point (i) follows. 

Point (ii) is quite delicate. First, there is C such that for all X e [1, oo), all y > 0, 
<^k{x)<Cx'' and - < C|x'' - 

Indeed, it is enough to prove that for Tk{x) — /q^(1 — cosG{z / x))dz, Tk{Q) = and |r^(a;)| < 
C. But Tk{x) = X J^^^'il - cosG{z))dz < x G^{z)dz, so that TKiO) = and |r'^(a;)| < 
J^{l-cosG{z))dz+x{K/x^){l~cosG{K/x)) < G^{z)dz + {K/x)G^{K/x), which is uniformly 
bounded by ([TTUl) . Consequently, for all X,Y eR^, 

\X<^>k{\X\) - Y^Ki\Y\)\ <C\X - Y\{\Xp + + Ci\X\ + \Y\)\\Xp - \Yp\. 
Using again that \x'^ — J/'*'! < 2|a; — y\/{x^^'^ + y^^^), we easily conclude that 

(3.2) \x<i>K{\x\) - y<^k{\y\)\ <c\x - + \Yn. 

Now we write 

= - [(u - -D) + (u* - u*)] • [{v - v^)^k{\v - - {v - w*)$Ar(|u - w*])] 
+ [{v -v) + (z* - w*)] • [(w - z*)$k(|w - z*|) - (u - - 5*1)] 

(3.3) =-[{v-i,) + (u* - u*)] • [{v - w,)$k(|w - w*|) - (w - z,)<I>k(|w - z.|)] 

+ (z, - u*) • [(w - z,)$if (|u - z*|) - (w - i;,)$ii:(|w - -5*1)] . 

By p.2p and the Young inequality, we deduce that 

Af <C{\v — v\ + |u* — u*|)|w* — z,|(|w — Uh.P + |w — z+P) 

+ C|z* — u*|(|u ~ v\ + |z, — w»|)(|w — z*!'"' + |w — u*P) 
<C[(|w — u| + |w, — v*!)^ + |w* — z,p(|t; — + \v ~ z^\^)'^] 

+ C\z^, — v^\{\v — -51 + |z, — u,| + — w*|)(|w — z^P + (|w — u| + |w — w,| + |w, — w*])'''). 
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The first term is clearly bounded by C{\v — -Dp + \v^ - + \vt, — + \v\ + + l^^*!)^''') 

which fits the statement, since 27 < 27/(1 — 7). We next bound the second term by 

C|z, — — z,| + |w — 

+ C|z, ~ ~ v\ + |w* — *^*|)''' 

+ C|z, — — v\ + |w, — — z*| + |w — u*!)'*' 

+ C|z, - - w| + |v* - ■5,1)-^+'''. 

Using that x^y'^ < x'^l^'^~^'> + (for the second line), that xyz'^ < (xz'^)'^ + (for the third line) 
and that xy^'^'^ < x'^/'^^~^^ + J/^, we obtain the upper-bound 

C|z, - w,p(l + \v\ + + |w,|)^ 

+ C{\v — v\ + |u, — u*!)^ + 1^* — w*P(|w — z^ \ + \v — w*!)^'' 
+ C(|w - w| + Iw, - v,\f + |z, - 
which is bounded by 

C{\v - 5p + |w* - + C|z, - w,p{(l + \v\ + \z,\ + \v,\y + \z, - w,|''/(2-T)-2 

+ {\V - Z,\ + \V- V,\f^ + \Z, - 

One easily concludes, using that max{7, 4/(2 — 7) — 2, 27, 2/(1 — 7) — 2} = 27/(1 — 7). 

We finally check point (iii). Using (|1.10l) . we deduce that 1 — cos(G(z/x''')) < G'^{z/x'') < 
Ciz/x^)-^/", whence ^'^(a;) < Cx^-y/" z'^/^dz = Cx^^/" K^-^/"" . Thus 

(3.4) Af{v, v^,v, C,) <C{\v - v^\^ + \v - v^\\v ~ i,4)\v - v^]^'''/'' K^-^/\ 

from which we easily conclude, using that \v ~ Vf,\\v — v^.]'^^^'^^'' <|u — p + — p"'"'''*'/''. □ 

We conclude with the hard spheres case. 

Lemma 3.4. Assume (11.31) . (|1.5p with 7 = 1, p.6p and adopt the notation of Lemma lS. 1\ 

(i) For all q > 0, there is Cq > such that for all M >1, all K £ [1,cxd), all G R'^, 

A^{v,v,,v,v,) <Af(|«-5p + |i;, -S,n + C,if(|^| + |«,|)e-*^'e^'(l^l'+l"*l'). 

(ii) For all q > 0, there is Cq > such that for all M > 1, all K £ [1, 00), all v, v-^,, v, G M^, 
(w, w*, u, u,) — A2 (v, z^,,v,Vt,) <AI{\v — -Dp + |w, — "5*1^) + C|w, — z,p(l + \v\ + \v^,\ + jz*])^ 

+ Cq{l + + |S,|)i^e-^'^'e'^'(l"l'+l"-l'+l^*l') 



(iii) For all q > 0, there is Cq > such that for all K Cz [1, 00), all v, w,, v, £ 
A^iv,v.,v,v.) <Cq{l + |i)|)e-^'e^'(l''l'+l-l'+l-l'). 
Proof. On the one hand, (jl.lip implies 



p3 



Ai {v,v^,,v,v^) <2ci\v — — 



\v - wJ)^ 



V ~ Vr \ + \v — Vt,\ 

<4c4(|w — u»| A |u ~ v*\){\v — + |w* — w*!^). 
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On the other hand, since G takes values in (0, 7r/2), we obviously have 
Consequently, we may write 

< M{\V - + \V.^ - V.^f) + Y^\^ ~ " '^*|l{4c4(|t)-l>.|A|-D-^.|)>Af}- 

Point (i) easily follows, using that |t; - v*\Ti^{4c4{\v~v,\A\ij~ij,\)>M} < \v - w*|1{4c4|^;-«.|>m} < 

Using all the computations of the proof of Lemma |3. 31 - (ii) except the one that makes appear the 
power 2/(1 — 7), we see that for Af^ :— (w, w*, w, v^,) — A2 {v, z*, -D, ?)*) 

Af <C[\v - vf + |u* - u*p + |v, - z^f{l + \v\ + + |z*|)^ + |z, - v*\{\v ~ + |w, - v^f)] 

<C(1 + |z» - - -Dp + |w, - w^p) +C|u, - 2;,p(l + \v\ + \v^\ + 

On the other hand, starting from p.Sp and using that 4>k{x) < t:K, we realize that 

Af < CK{1 + \i\ + + Ivp + |w,p + |z,p). 

Hence we can write, for any M > 1, 

Af <M(|v - up + - u,p) +C|u, - z,p(l + |t;| + |u,| + |z,|)2 

|)(l + |up + |u*P + |z,P)l 

But (1 + |wp + |w,p + |z*P)a{C(l + |z.-i).|)>Af} < (1 + |v| + + \z^\)'^l{c{l+\v\ + \v,\ + \z,\)>M} < 

(1 + \v\ + \v,.\ + |2,|)2e-^/'eC-(i+|H+k.|+|z.|)' < ^^g-A/'gC,(|.r+|..r+|..r)_ poj^^- (jj) checked. 

Finally, we observe that "^^ix) < tt G^{z/x)dz. But here, G'(z) = (^/2 — z)+ whence 
'I'i<-(a;) < {tt'^/24)x1[x>2K/7t} < 5.Ta{j.>x/2}- Thus for any g > 0, ^^(a;) < 5a;e^^''e^'''^', so that 

A^{v,v,,i,i,) <C{1 + |5|)(1 + Ivp + |t;,p)e-^> - w^le^'l----!" 

as desired. □ 

4. Convergence of the particle system with cutoff 

To build a suitable coupling between the particle system and the solution to we need 

to introduce the (stochastic) paths associated to To do so, we follow the ideas of Tanaka 

[551 131] and make use of two probability spaces. The main one is an abstract (51, Pr), on which 
the random objects are defined when nothing is precised. But we will also need an auxiliary 
one, [0, 1] endowed with its Borel a-field and its Lebesgue measure. In order to avoid confusion, 
a random variable defined on this latter probability space will be called an a-random variable, 
expectation on [0, 1] will be denoted by E^, etc. 

4.1. A SDE for the Boltzmann equation. First, we recall the classical probabilistic interpre- 
tation of the Boltzmann equation initiated by Tanaka [351 [SH] in the Maxwell molecules case. 

Proposition 4.1. Assume ((TH) or and let fo £ P2(R^)- Ifl 6 (0, 1], assume 

additionally (jl.8l) and that /q satisfies (|1.12p . Let {ft)t>a be the corresponding unique weak solution 
to (II. ip . Consider any /q- distributed random variable Wq and any independent Poisson measure 
M{ds,da,dz,d(p) on [0,oo) x [0,1] x [0, 00) x [0, 27r) with intensity measure dsdadzdip. Consider 
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also, for each t >0, a ft- distributed a-random variable , in such a way that (i,a) Wf{a) is 
measurable. Then there is a unique (cddlag adapted) strong solution to 



(4.1) Wt=Wo+ / / / c{Ws-,W:{a),z,^)M{ds,da,dz,d^). 

Jo JR3 Jo Jo 

Furthermore, Wt is ft- distributed for each t>0. 

We will note {Wt)t>o such a Boltzmann process. It can be viewed as the time-evolution of the 
velocity of a typical particle in the gas. 

Proof. The proof is very similar to that of fTT, Proposition 5.1], see also [TH Section 4] and 
is omitted. In 117', Proposition 5.1], the same Boltzmann equation is studied, with much less 
assumptions on /g (so that uniqueness is not known for (II. ip ). But the formulation of the SDE 
is different (it is equivalent in law). The same proof as in jl7] Proposition 5.1] works here, with 
several difficulties avoided due to the facts that /o has exponential moments and that uniqueness 
is known to hold for (II. ip . □ 

4.2. A SDE for the particle system. Here we write down a Poisson stochastic differential 
equation corresponding to Nanbu's particle system and we prove Proposition 1 1 . 21 - (i V 

Proposition 4.2. Assume dO]), dm]), dHH) or dUT]) and let fo e P2(R^), N > 1 and K £ 
[l,oo). Consider a family {Vf!;)i=i^,,,N of i.i.d. f^- distributed random variables and an independent 
family {Of {ds,dj,dz,d(p))i=i^,,,^N of Poisson measures on [0,oo) x {1,...,7V} x [0, oo) x [0,27r) 
with intensity measures ds (^N~^J2^=i^k{dj)j dzdif. There exists a unique (cddldg and adapted) 
strong solution to 

(4.2) ^/■^■^ = T/J+/ // / CK(V;^''',F/::^''',z,^)Of(ds,dj,dz,d^), i^l,...,N. 



Jo JjJO JO 

Furthermore, {Vt^'^'^)i=i^,,,^N,t>o is Markov with generator Cn,k- We have 



/r3 H'^foidv) and, if J^^ \v\Pfo{dv) for some p>2, sup[o,T] ^ 



< a 



P.TJo,K- 



Proof. First of all, observe that we actually deal with finite Poisson measures, since ck vanishes for 
z > K. Thus, strong existence and uniqueness for (14. 2p is trivial: it suffices to work recursively on 
the instants of jumps (which are discrete) of the family {Of^ {ds,dj,dz,d(p))i=i^,,,^N- Consequently, 
Vf = {Vt^'^'^ , . . . , Vj^'^'^) is a Markov process, since it solves a well-posed time-homogeneous 
SDE. Its infinitesimal generator is classically defined by (|2.7p . with actually a sum over all cou- 
ples £ {1, . . . , N}^ , but this changes nothing since the terms with i ~ j vanish because 
ck(v,v, z,tp) = for all u e R^. Next, a simple computation shows that 



1 A^ poo p2Ti 

nvt'^^^^?]=nv,'?] + ^,yi / / ^{\vt''-^+CK{vy-'',vi-^-^,z,^)\^ 



Nj-,Jo JO 



\Vy'''f)d^dzds 



N - 1 



=E[|y„in-K^ / / / E{\civy'^,vy'^,z,^)\' 







+ 2Vl'«'^ . ciVy^'<, ^)) d^dzds 
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by exchangeability. But, as seen in the proof of Lemma |3.1[ 

[ \ c(v, w*, z, + 2t; • c(w, w*, (p)]dipdz — [|w — — 2v ■ {v — — 



whence, using again exchangeabihty. 



l,Af,K|2i 



=E[|K,T] 



=n\Vo'\'] 



N ~1 
N 



N -1 



l,N,K 



l.N.K 



2.N.K\ 



1,N.K _ y2,N,K 



\)\ds 



N 



yl,N,K _ y^ 
2,N,K ni2,N,K 



2,N,K\2 



l.N.K 



1,N,K 



^yl^M^K 
1,N,K 



y2.N.K^ 



In this last expression, the integrand is zero, so that, as claimed, E[|V/'™'"p] = E[|Vj3^p] = 



^,N,K^2^ 

/r3 \v\^fo{dv). Recalling finally (I2.ip and (|2.5p . we see that |c(w, w,, z, (y9)| < |t; — i;, 
P> 2, 

(|w + c(«, z, - \v\P)d^dz < CpK{\v\ + 



Thus for 



Consequently, we obtain as previously 



and conclude, using again exchangeability, that E[|y/'^'-^|P] < E[|y(/|^']e^'^p-'^* as desired. □ 

This allows us to deduce 

Proof of Proposition The strong existence and uniqueness for the SDE (|4.2p classically 

implies the existence and uniqueness of a Markov process with generator Cn,k- D 

4.3. The coupling. Here we explain how we couple our particle system with a family of i.i.d. 
Boltzmann processes. For example, we want to couple V^''^''^ with a Boltzmann process W^. The 
main difficulty is that at each collision, is collided by an independent particle (using W^) while 



' ' is collided by some V/' ' . We thus have to choose j in such a way that V^' ' is as close 
as possible to W^, but j has to remain uniformly chosen. 

A technical problem obliges us to introduce the set (M"^)^ := {w G (R^)''^ : Wi wj \/ i ^ j}. 

Lemma 4.3. Lei /f e C([0, oo), 7^2(1^^)) be such that ft has a density for all t > . Let also N > 
1 be fixed. For v = (vi, . . . ,vn) & 



N 



we denote by fi^ := N Sy- the empirical 



associated to v. There exists a measurable map {t, w, v, a) i-^ (VFj*(a), Z^(w, a), Vf*{v, w, a)) from 
(0, oo) X (R'^)^ X (R'^)^ X [0, 1] into R'^ x R'^ x R^ enjoying the following properties 

(a) for all t > 0, the a-law of is ft, 

(b) for allt>0,w € (R^)^, the a-law of Z;(w, .) is fi^, 

(c) for allt>0,w e (R^)^ , v G (R^)^, the a-law ofVt*{v,vf, .) is , 



(d) for allt>{),we 



N 



V G 



N 



, thea-lawof{Z;{w,.),V;{v,w,.)) is N ^ Ef '^(«..,«.)' 



(e) for all t>0,all^ve (R^)^, \Wt*{a) - Z;{^v,a)\^da = 
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Proof. We first consider, for each t > 0, Wf such that point (a) holds true and such that {t, a) i— )■ 
W^{a) is measurable. 

Next, we recall that by Brenier's theorem (see e.g. Villani [42, Theorem 2.12 p 66]) for each 
t > and each w G (M'^)^, since ft does does not charge small sets (because it has a density 
by [H]), there exists a unique map i^t,w : M'^ n> R'^ such that, setting Z^{w,a) := Ft_w(Wj*(a)), 
points (b) and (e) hold true. In other words, (M^f*(.), Zj*(w, .)) is an optimal coupling for ft and 
/i^. Furthermore, Fontbona-Guerin-Meleard [16^ have shown that Ft^-w{x) is a measurable function 
of (t,w,x). Consequently, Zj*(w, a) is a measurable function of (t,w,a). 

Finally, we define, for any w G (K'^)^ and any v G (R'^)^, the map Gw.v '■ {wi, . . . ,wn} h> 
{vi, . . . ,vn} by Gw,v(u;i) = Vi (here we need that w G (M^)^). We then we put Vt* {•v,-w,a) — 
Gw,v('^t*(w, a)), which is clearly measurable (in all its variables). Point (d) follows from (b) and 
the definition of Gw,v and finally (c) follows from (d). □ 

Here is the coupling we propose. 

Lemma 4.4. Assume ((Ot . (fTSl) . (fTG)) or (fTTI . Let fo G 7'2(R^)- Assume additionally (fTSl 
and ()1.12p if ^ ^ (0, 1]. Lei {ft)t>o be the unique weak solution to and assume that ft has a 

density for all t > (see Theorem ] 1.1]) . Consider N > 1 and K G [l,oo) fixed. Let (y^^)i=i^...N 
he i.i.d. with common law fo and let {Mi{ds,da,dz,d(p))i=i^....N be an i.i.d. family of Poisson 
measures on [0, oo) x [0, 1] x [0, oo) x [0, 27r) with intensity measures dsdadzdip , independent of 

iVi^)i=l....N- 

(i) The following SDE's, for i = 1, . . . , N , define N independent copies of the Boltzmann process: 

Wi=V^+ / / / c{W:_,W:{a),z,^)M,{ds,da,dz,d^). 
Jo Jo Jq Jo 

Ln particular, for each t > 0, {Wt)i=i^...^N o.fe i.i.d. with common law ft. Consequently, since ft 
has a density for all t>Q, {Wl)i=i^,,,^N G (R^)^ a.s. 

(ii) Next, we consider the system of SDE's, for i = 1, . . . , N , 

^/'^■^ = 1^'+/ / / / CK(F;i^^^,l/;(Vf_:^,W,_,a),z,^ + ^,,„,,)M,(ds,da,dz,d^), 
Jo Jo Jo Jo 

where we used the notation Vf_:^ = (F/i^'^, . . . , Vl['^'^) G (R^)^, W,_ = {W}_, . . . ,Wfi) e 

(R3)^ and where we have set ip^,a^s ¥'o(W^l- - W:{s,a),V,'L^''^ - V;{V^1^ ,Ws-,a)) for 
simplicity. This system of SDEs has a unique solution, and this solution is a Markov process with 
generator Cn,k md initial condition (Vo*)i=i^...jv- 

(Hi) The family {{Wt ,Vl'^'^)t>o, {W^ ,Vt^'^'^)t>o) is exchangeable. 

Proof. Point (i) is a direct consequence of Proposition l4. 1 1 and point (iii) follows from the exchange- 
ability of the family (Vq*, Mi)i=i^...^jv and from uniqueness (in law). In point (ii), the existence and 
uniqueness result is also immediate, since the Poisson measures under consideration are finite (or 
rather, are finite when z is restricted to [0, i^], which is the case since ck = cTL^^kk})- Finally 
{Vt^'^'^ , . . . ,Vt^'^'^)t>o is a Markov process with generator Cn^k due to the fact that for all 
V G (R^)^, all w G (R^)^, all s > 0, all (pij G [0, 27r), for all bounded measurable function 
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E 



i=l 
N 



JO JO 



-E 



OO /'27T 



(^(f){-v + CKivi,V*{-v,w,a), z,ip + ipij).ei) - 4){-v)^ dipdzda 
{4){-v + Cif(ui, , z, <p + Lpij).ei) - 4)(y)j dipdz 
(^(f){-v + CKivi,Vj,z,ip).ei) - 4)(y)j dipdz, 



which is nothing but Cn,k4'{^)^ see (j2.7p . We used Lemma [4.3K c) for the first equahty and the 
27r-periodicity of ck (in ip) and the fact that CK{vi, Vi, z,p) ^ for the second one. □ 

4.4. Estimate of the Wasserstein distance. We can now prove our main result in the case 
with cutoff. We first study hard potentials. 

Proof of Theorem\r^(ii) when K e [I, oo). Wc thus assume (fO)) . ([T3|) with 7 G (0, 1) and ([TTT)) . 
We consider /q £ 7'2(K'^) satisfying (|1.12l) for some p G (7, 2) and fix g G (7,p) for the rest of the 
proof. We also assume that /o is not a Dirac mass, so that ft has a density for alH > 0. We fix 
> 1 and K £ [1, 00) and consider the processes introduced in Lemma [4.41 

Step 1. A direct application of the Ito calculus for jump processes shows that 



(^t /•! /"OO /'27r 

E 



Jo Jo Jo 



\Wl - V}'^'^ + l\'^{s,a,z,ip)\'^ - \Wl ~V}'^'^\'^]dipdzdads, 



where 



Using Lemma |3.1[ we thus obtain 

JO 

where, for i = 1, 2, 3, 



da. 



5tep 2. Using Lemma [3. 31 - (i). we see that for all M > 1 (recall that q G (7,^) is fixed) 

-Bf (s) <M (\[\Wl - F/'^'^p + \W:{a) - V;*(Vf'^, W„a))p] da 
Jo 

+ Ce-*^''" / E [exp(C(|VF,y + |VF;(a)|'?))] da 
"'0 

<M [\[\Wl ~ F/'^'^l^ + \w:[a) - K*(V^^,W„a))l2] da + Ce"^^"' 

JO 
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To get the last inequality, we used that W} and W*(.) are independent and satisfy ^ fs and 
W*{.) fs, whence 

[\[exp{C{\Wl\'> + \W:{a)n]da = ([ e^l™' Vs(du;)) ' < oo 

Jo VH3 / 

by (ma. 

Step 3. Roughly speaking, B2 should not be far to be zero for symmetry reasons. We claim 
that B2 would be zero if W* [a) was replaced by Z* (Ws , a) . More precisely, we check here that 

B^{s) E [Af (T^i, Z:(W„ a), F/'^'^, ^/(Vf W„ a))] da = 0. 



By Lemma l4.3l -(d). we simply have 



i,N,K-. 



i=l 



N - 1 
N 



2,N,K\ 



by exchangeability and since A2 (w, w, w, w) = 0. Finally, we write, using again exchangeability. 

This is zero by symmetry of : it holds that A^{v, w,, v, w,) + A2 {v, w*, v, u*) = 0. 
Step 4- By Step 3, wc thus have 



E [Af{W^, W:{a), Vy^"", V;*(Vf ^, W„ a)) 

- A^iWl Z:(W„ a), P;!-^'^, F;(V^^, W„ a)) 
Consequently, Lemma l3.3l -(ii) implies 



da. 



B^{s) <C / E \W^ - y/'^^-^l^ + \W:{a) ~ l/;(Vf W„a)|^ 
Jo 

+ \W:{a) ~ Z:(W„a)p(l + \Wl\ + \W:{a)\ + |Z:(W„ a)|)2^/(i-^) 
Step 5. Finally, we use Lemma 13.31 - (iii) to obtain 



da. 



da. 



Since W} ~ J^, we deduce from (11131) that E[|M^,1|4t/''+2] ^ jj^^ |^|47/>^+2 < By Lemma 
we also have W;{.) - whence /^^ |VK;(a)|''^/'^+2da = /j^a < C. Propo- 

sition 112] shows that E[|t//'^'^|2] = /jj3 |i)|2/Q(di,). We next infer from Lemma Olfc) that 
/o lK*(Vf^,W„a))|2da = iV-i^f IK^'^'^'I'- Consequently, E[/„^ |K*(Vf W„ a))|2da] = 
My}'^'^?] = /r3 |f P/o(«^t')- As a conclusion. 



(s) <CXi-2/^ 
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Step 6. We set u^'^ := IE[| M^/ — K^'^'^P]- Using the previous steps, we see that for all M > 1, 



N,K 



<Cte 



+ CtK^-^/'' + {M + C) f [uf'^+f E[\W:{a)-V;iV^''\Ws,a)\^]da]ds 
Jo Jq 



C 



E 



\W:{a) - Z:(W„a)p(l + \wn + \W:{a)\ + |Z:(W„a)|)2^/(i-^) 



dads. 



We now write, using Minkowski's inequality and Lemma l4.3l -(d) and (e), 

1/2 



(4.3) 



|iy;(a)-y;(vf^\w,,a)| 



da 



< 



\W:{a)-Z:{W,,a)\ 



da 



1/2 



E 



1 ^ 



|z;(w„a)-F;(v^^,w„a)p 

1/2 



da 



1/2 



=E[W|(/.,/z^J]i/^ 

=E[W|(/.,M^J]V2 + (^.f^)V2 
by exchangeability. We deduce that 

/ ' E[\W: {a) - V; (Vf W, ,a)\']da< 2E[W| (/, , fi^^ )] + 2u 



(4.4) 



Next, a simple computation shows that for all e G (0, 1), 



(4.5) 



E 



< 



< 



\W:{a) - Z;(W«,a)|2(l + \W}\ + \W:ia)\ + |Z:(W„a)|)T^ 
\W:{a) - Z:(W„«)|2-(1 + \Wj\ + \W:ia)\ + \Z:(W,,a)\)-- 



da 



da 



\W:ia)-Z:iW,,a)\' 



E 



da 



(1 + \W}\ + \W:{a)\ + \Z;{Ws,a)\)^o^^^ 

<C,(E[W|(/„M^J])'^. 
For the last inequality, we used Lemma r4.3l -(e). the fact that by (I1.13p . 



da 



E 



\W}\ 



+2 



, . 4-2 



\W;{a)\^o^^^da= I \v\^o^+^ fsidv) < 



and that, by Lemma B751 -fb) 

/' 



E 



l^:(W.,a)| 



r+2 



da =E 



1 ^ 



I e(l-T) 



+2 



E 



+2 



< a 



22 



NICOLAS FOURNIER AND STEPHANE MISCHLER 



We end up with: for all e e (0, 1), ah M > 1, 

Jo 







Now we observe that E[yV'|(/s, /^^J] = SNift), recall pTTil) . because W} , W^/^ are i.i.d. and 
/(-distributed. Since £jv(/t) < 2 /j^a [^^/(((iu) = 2 /j^g |u|^/o(dw), since M > 1 and K G [1, oo), we 
get 



where we have set 



SN,K,t :=i^i-2/'' + supejv(/s). 

[o,t] 



Hence by Gronwall's lemma, 



sup «f <C,T fe-*^"'^ + MS]if^ e^'*^^, 

[O.T] ^ ' ' ' 

this holding for any value of A/ > 1. We easily conclude that 

N.K ^f-, rl-e 

[0,T] ■ ' 

by choosing M = 1 if 5i^^k,t > 1/e and M = | \og5N,K,TV^'^ otherwise, which gives 

sup T/f <C, (t6n,k,t + S],-'Jl\ \ogdN,K,T\''^') e^^l'°s«"'-'-l'''^ < C,,t4"^t> 

[0,T] ^ ' ' ^ 

the last inequality following from the fact that j/q < 1. 
Final step. We now recall that /if = I^-^n.k and write 

But E[W|(/i;;^«.., /i^J] < E[7V-i Ef l^/'"^'"" " W^^] = E[|y/'^'^ - W'.ip] = u^-"" by exchange- 

ability, and we have already seen that E[VV|(/t^^ , /t)] — e7v(/t)- Consequently, for all e G (0, 1), 
aUie[0,T], 

E[w2(A^f < a,T<5]v"^T + 2eiv(/0 <Ce,T (k^-^/'' + sup eNift)] 

and this proves (|1.18p . Using finally (jl.lSp and applying finally Theorem 11.31 (with any choice of 
fc > 2), (fTTO)) easily follows. □ 

We next study the case of Maxwell molecules. 

Proof of Theorem \1.4\ -(i) when K g [l,oo). We thus assume (|1.3p . (jl.Sp with 7 = and (|1.7p . 
We consider /o G T'2{^'^) not being a Dirac mass. We also assume that /o G P4(M^) or that 
/r3 /o(^^) log/o(w)'iw < 00, so that ft has a density for all t > 0. We fix > 1 and iiT G [1, 00) and 
consider the processes introduced in Lemma 14.41 
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Step 1. Exactly as in the case of hard potentials, we find that 



l.N,K\2- 



< [ [B^{s)+B^{s) + B^{s)]ds, 
Jo 



where Bf{s) := />[Af (PF^ T^;(a), F/^^^^, K*(V^^, W„ a)) 
Step 2. By Lemma E^Ki) , we have 5f (s) = 0. 



da for 1 = 1,2, 3. 



Steps 3 and 4- By Lemma [23]- (ii), it holds that for = /o^(1 ~ cosG(z))d. 



da. 



Jo 

Step 5. By Lemma IX^ fiii) 

(s) < CK^-^^" [\\\W^\^ + \W:{a)\^ + \Vy'^\^]da < CK^-^'" , 



since, as usual, E[|PF,i|2] = \W*{a)\^da = E[|y/'^^-^|2] = j^^ \v\^fo{dv). 
Step 6. Setting uf := £[1^^^! - T^^'^'^l^], we thus have 



N,K 



<CK^-^/''t + C 



K 



,N.K 



|w^;(«)-v;*(vf^,w„a) 



da ds 



by (|4.3p . Next we recall that eM{ft) — E[W|(/s, /^wj]' we set Eat^t = sup[o,T] ^nUi) and we recall 
that (k < /o°°(l ~ cosG{z))dz < oo. We thus may write, for aU t G [0,T], 



Jo 



N,K 



Then we have (wf'^)' < Celj^^iv^-'^y/^ , so that (wf''^)!/^ < (C(i^i-2/.^g^_^j.)y)i/2^^£V2^^_ 
We conclude that 

sup wf '-^ < C(ii:^"2/'^ + EN tT)T + CT^ea, t < CK^-^^^T + C{T + T^)eN t- 

[0,T] 

Final step. Exactly as in the case of hard potentials, for t G [0,T], 

E[>V|(/if /*)] < 2eNift) + 2uf < CK'-'/'T + C{1 + Tf sup £jv(/0 

[0,T] 

whence (I1.16p . If finally /o G 7^^(111^) for aU fc > 2, then we know that sup[g /j^a \v\'^ ft{dv) < oo, 
so that (jl.l7p follows by application of Theorem 11.31 □ 

We conclude with hard spheres. 

Proof of Theorem \1.4[ (m). We thus assume (|1.3p . (|1.5p with 7 = 1 and (|1.6I) . We consider 
fo G 7^2 (R'^) satisfying (|1.12p for some p G (7, 2) and fix g G (7,p) for the rest of the proof. We also 
assume that /o has a density, so that ft has a density for all t > 0. We fix > 1 and K G [1, 00) 
and consider the processes introduced in Lemma 14.41 
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Step 1. Exactly as in the case of hard potentials, we find that 



u. 



N,K 



< f [B^is) + B^is) + B^is)]ds, 
Jo 



where Bf (s) := /> Uf W^;(a), F/'^^^, ^/(V^^, W„ a)) 



da for 1 = 1,2, 3. 



Steps 2, 3, 4, 5, 6. Following the case of hard potentials, using Lemma [3.41 instead of Lemma 
1 we deduce that for all M > 1, 

^i?f (s) <2M l\[\Wl ~ F/'^^^'P + \W:{a) - K*(V^^,W„a)p] da 



+ C E 



E 



(l + |Ki-^-^| + |K*(Vf'^,W„a)|) 



X e 



c( I VK,' r + 1 (") r + 1 2; ( w, ,q) I " ) 



da 



|w^;(a)-z:(v^^,w„«)|2 

X (1 + \W}\ + \W:{a)\ + |Z:(Vf ^, W„a)|)2' 
Proceeding as in (|4.5p . we deduce that the last line is bounded, for all e G (0, 1), by 

and using (I4.4p . the first term is bounded by 

4ME [W2(/„^^J] +6Muf'^. 



da 



Using finally the Cauchy-Schwarz inequality, that, thanks to Lemma [4.3l -(c) and by exchangeability, 
E[/o lK*(Vf'^,W„a)pda] = E[iV-i IK^^^^^P] - = J^, \v\' fo{dv) < oo and 



(|1.13p . we easily bound the second line by C(K: 



(recall that - /«, that W*{.) fs 



and that, by Lemma|13(b), e^l^= (^-"^I'da = N'^ J2i e^l^=l'). 

Recalling that E [W|(/s,^^J] = £jv(/s) and setting eAr^t ~ sup[o_t] eA'(/s), we thus have, for 



any M > 1, any e G (0, 1), 



N,K 



<6M [ u^'^ds + Ct{Ke 
Jo 



Thus by Gronwall's Lemma, 



N,K 



<Cet{Ke 



Choosing Af = 2/^ and using that Ke~^^^'>'' < Ce~^', we deduce that 

sup «f <C,T(e-^' + 4"T^')e^^^^ = C,T(e-^' + (sup eM{fs))'-'/')e''^^ ■ 

[0,T] ' [0,T] 



Final step. We conclude as usual, using that E[>V|(^f '^,/t)] < 2eN{ft) + 2wf'-^ to obtain 
([L20)) and then (|LT3)) and Theorem O to deduce □ 
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5. Extension to the particle system without cutoff 

It remains to check that the particle system without cutoff is wcll-posed and that we can pass 
to the Hmit as i^T — >• cx) in the convergence estimates (|1.16p - (fTTT7l) - (|1.18"l) - (|1.19p . We will need the 
following rough computations. 



Lemma 5.1. Assume (|1.3p . (|1.5p and (|1.6p or (|1.7p . Adopt the notation of Lemma \3.1\ There 
are C > 0, k > and (5 > (depending on ^,v) such that for all K G [l,oo), all G M'^, 

3 

^Af{v,v^,v,v^) < C(l + + + \i\ + \v^\Y{\v - + - + K^^). 



1=1 



Proof. Concerning A^ , we start from p.ip (this is valid for all 7 G [0, 1]) and we deduce that 

(u, u*, u, u*) <8c4(|u — -Dj A |u, — v,,\Y{\v — v\ + |u, — w*])^ 



<C{l + \v\ + \v,\ + \v\ + \i,\f{] 



We then make use of (|3.2p (also valid for all 7 G [0, 1]) to write 

A2 {v, u*, w, w*) <C{\v — v\ + |w* — 5*1)^(1?; — w*p + |w — 5, p) 

<C(1 + |w| + + |w| + |w*|)^(|w - + |w* - 

For , we separate two cases. Under hypothesis (|1.7p . we immediately deduce from p.4p that 

(i;,v,,C,C,) <C(1 + |w| + \v,\ + \i\ + |f},|)2+2T/-i^i-2/.^ 

Under hypothesis (|1.6p . we have seen (when 7 = 1, at the end of the proof of Lemma [3.4p that 
'^k{x) < 5a;^ll{^7>;<-/2}j whence k{x) < IQx^'^ /K and thus 

Af (w,u,,w,w,) <C{\v - v,\y \v - i,\f+'^^ K-^ < C{l + \v\ + \v,\ + \i>\ + \v^\)^+^m~\ 

The conclusion follows, choosing k = 2 + and 6 = 2/v — 1 under (|1.7p and k = 2 + 27 and 
(5 = 1 under (dH). □ 

Now we can give the 

Proof of Provosition \1.2\ -fii). We only sketch the proof, since it is quite standard. In the whole 
proof, > 2 is fixed, as well as /o G 7^2 (K'^) and a family of i.i.d. /o-distributed random variables 

{Vq )i=l,...,N- 

Step 1. Recall (12.81) . Classically, (l^*'^'°°)i=i,...,Ar,t>o is a Markov process with generator Cn 
starting from (V^*' )i=i N if it solves 

ft /> /"Oo />27r 

(5.1) = K + / // / c(t/;:^^°°,Kf::^^°°,z,^)of(ds,dj,dz,dv.), z = i,...,iv 

Jo Jo Jo 

for some i.i.d. Poisson measures [ds, dj, dz, d(/?))i=i,...,Ar on [0, 00) x {1, . . . , N} x [0, 00) x [0, 2n) 
with intensity measures ds (^N^^ J2^=i ^kidj)j dzdtp. 

Step 2. The existence of a solution (in law) to (|5.ip is easily checked, using martingale problems 
methods (tightness and consistency), by passing to the limit in (|4.2p . The main estimates to be 
used are that, uniformly in K ^ [1, oo) (and in > 1 but this is not the point here). 



E[|F/^^^^n = / \v\'f^{dv) and E 

R3 



sup \V^ ' ' I 

[0,T] 



< Ct 
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for all T > 0. This second estimate is immediately deduced from the first one and the fact that 
Io° lo^ |c(z;, w*, z, (/?)| < C\v — v*\^~^'^ < C(l + \v\ + \v^\)^. The tightness is easily checked by using 
Aldous's criterion [T]. 

Step 3. Uniqueness (in law) for (|5.ip is more difficult. Consider a (cadlag and adapted) solution 
(V/' '°^)i=i....^N,t>o to (|5.1I) . For K e [l,oo), consider the solution to 

T//'^'^ = ^'+/ // / CKiV:L''''',V,t''''',z,v + ^s,^,J)0^{ds,dJ,dz,d^), z = l,...,iV 
"'0 Jj Jo Jo 

where ips,i,j ■= 'Po{Vs-^''°° — vjl^'°° ,V^!:^'''^ — V^!^'^). Such a solution obviously exists and is 
unique, because the involved Poisson measures are finite (recall that ck{v, u*, z, 1^9) = for z > K). 
Furthermore, this solution (V/' ' )i=i^...^N,t>o is a Markov process with generator Cjq.K starting 
from (Vj'^)i=i_..._Ar (because the only difference with (|4.2I) is the presence of ^Ps.i,j which does not 
change the law of the particle system, see Lemma l4.4l -(ii) for a similar claim). Hence Proposition 



11.21 -fi) implies that the law of (V/' ' )i=i....,Ar,(>o is uniquely determined. 

We next introduce tn,k,a = inf{i > : 3 i e {1, . . . , A^}, \V^-'^-''^\ + |V/'^'-^| > A}. Using, on 
the one hand, the fact that (V/'^'°°)i=i^...^7v,t>o is a.s. cadlag (and thus locally bounded) and, on 
the other hand, the (uniform in K) estimate established in Step 2, one easily gets convinced that 

(5.2) V T > 0, lim sup Pr[rAr ka<T]^Q. 

K>1 

Next, a simple computation shows that 



Trl,N,K |2i ^ ^ Y^Tirr / / / Mt/I.A'.oo ^rl.N.K , ^l.j,N,K , 

V^tXr™,^,^! ] <T^2^IE[ / / / [\V,_ ^ -V,_ +A,^' ■ (z,^)| 



where 

f.l,j.N,K, N /yrl.Af.OO -.rl.N.OO \ /^rl.N.K jrl,N,K , N 

Using Lemmas 13.11 and 15.11 and the fact that all the velocities are bounded by A until T]y x,A, we 
easily deduce that 

jgn^l^AT^OO _yl,N,K |2l 

N Uo 

<Ct{1 + AfK-' + C{1 + AT f , - F/xf/, , ^ds 

by exchangeability. We now use the Gronwall lemma and then deduce that for any A > 0, 

lim I 

Gathering (15.21) and (|5.3p . we easily conclude that for all t > 0, v/'^'^ tends in probability to 
pri.w.oo as AT —j> 00. Thus for any finite family < < ••■ < t;, {V^".^"^)i^i^,,,^N,j=i,...,i goes 
in probability to {V^^^'^)i=i^,,,^Nj=i^,,,^i, of which the law is thus uniquely determined. This is 
classically sufficient to characterize the whole law of the process (V/'^'°°)i=i_...^Ar_t>o. 



(5.3) jim_ sup E[|FAv^-,, - F/xf;^,, p] = 0. 
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Conclusion. We thus have the existence of a unique Markov process (V^''^'°°)i=i,...,Ar,t>o 

with generator Cn starting from {VQ'^)i=i n, and it holds that for each t > 0, each N > 2, 

{V^'^'°°)i=i^,,,^N is the limit in law, as iiT ^ oo, of {Vf'^'^)i=i^,,,^N- O 

To conclude, we will need the following lemma. 

Lemma 5.2. Let N > 2 be fixed. Let (^''^'^)i=i....,jv be o, sequence of (R.^)^ -valued random 
variable going in law, as K ^ oo, to some {E?)^ -valued random variable (-f '^)i=i,...,jv- Consider 
the associated empirical measures v^'^ :— '^^^■^6xi.N,K and :— N^^^^^^Sxi.^ . Then 
for any <? S p2(IR^), 

E [Wl < liminf E [>V| • 

Proof. First observe that the map {xi, . . . , xn) m> W2{N^^ 6xi, g) is continuous on (R'^)^. 
Indeed, it sufhces to use the triangular inequality for W2 and the easy estimate 

/ N N \ N 

Consequently, W| (^i/^'^,g) goes in law to W| (^v^ ,g). Thus for any A > 1, we have 

E [Wl {i^^,g) A A] = ^hrn^E [Wl A A] < liminf E [>V| • 

It then suffices to let A increase to infinity and to use the monotonic convergence theorem. □ 
This allows us to conclude the proof of our main results. 

Proof of Theorem \r^(i)- (a) when K ^ 00. Recah that (fni| - (fTTf| - (fn^ - (frT9l) have already 
been established when K G [l,oo). Since (V^*'^'"")^^!,...,^ is the limit (in law) of (V/'^'^)i=i^...^jv 
as ivT 00 for each t > and each N > 2 (see the conclusion of the proof of Proposition II .21 - fii)) . 
we can let if 00 in ([TTT6| - ([TTf| - ([TT^ - (frT9l) using Lemma [O □ 

Appendix A. Quantitative empirical law of large number 
We conclude the paper with the 

Proof of Theorem \1.3[ We thus assume that / G Vki^'^) for some k > 2, with d — 3. We denote 
by fi^ the empirical measure associated to N i.i.d. /-distributed random variables. We denote by 
Mkif) := J^, \v\''fidv) and observe that E{Mk{p^)) = Mk{f). 

Step 1: k > 20. Consider the Sobolev norm defined, for /, 5 G ViW^), by 
ll/-.9lll^-. / {l + \i?)-'\f{^)-m)?di, 
where / and g denote the Fourier transform of / and g. We have (see [31] Lemma 4.2]) 

E(||m^ - /111,-.) + - l/(?)l'K < %^ if ^ > d/2, 

as well as the classical inequality (see e.g. [25l Lemma 2.1]) 

W2{f,g)<CdAMkif + 9))^\\f-g\\^ if s>d/2. 
In dimension d — 5 with the choice s = 3/2 (which is forbidden), this gives 

mif^g) < C{Mk{f + g))^\\f ~ g\\f0!,, if s > d/2. 



28 



NICOLAS FOURNIER AND STEPHANE MISCHLER 



Consequently, by Holder's inequality, 



2k-4 

-3 
3/: 



k~2 
fc + 1 



<CMfc(/)^E(|!A.^-/||^_3/.)^ 
<CMfc(/)^iV"^. 

Since the choice s = (i/2 is forbidden, we only get an estimate in N^"^"^ for any r < (fc — 2)/ (6fc + 6). 

Step 2: k € (7, 20]. We sketch the proof of Rachev-Ruschcndorf result [35l Theorem 10.2.1] to 
which we refer for details and that we will slightly modify in the next step. Let 

/.^■^ := * r ■■= f * .9^ fix) e-^g{x/e), g{x) = {2^)-^''' e-\^\"l\ 
Using a standard L^-upperbound of W| and the Cauchy-Schwarz inequality, 

<2de' + 2\\\x\'{f,^'%x)^r{x))\\L. 

<2de^ + 2||(1 + |a;|)2+('^+")/2(i + \x\y^^+''^/^{^i^'%x) - fix))\\Li 
<2de' + C„||(l + \x\)^+''+^{fi''''{x) - r{x)f\\]^,' 
for any a > 0. By Jensen's inequality, we get, for any e > 0, any q > 0, 

/ \ 1/2 

< 2de' + C„e(||(1 + \x\r+^+^{fi^^%x) - r{x))'U.) . 
But this last quadratic quantity can be easily computed and one gets 

E (11(1 + \x\)^+^+^{fi^^-{x) - nx))'\\L^) - ^ / (1 + * f){x)dx. 

Since d — S and since / e VkiM.^) with k > 7 ^ 4 + d, we may choose a = k — 7 and we easily 
deduce that for all e e (0, 1), /^.(l + \x\)^+'^+°' (g^^ * /)(a;)dx < Ce-3(i + Mfe(/)). We thus have 



E(W2 (a^'\/)) <6e' + C- 



iVl/2£3/2 

for all £ e (0, 1). Choosing e = iV-^/^, we find E(W|(Ai^, /)) < CN-'^''^ as desired. 
5tep 5; k G (2, 7). We start as in the previous step and write 

Wlin^J) <2de' + 2\\\x\\^,''''{x) - r{x))h. 

<2de' + 2||(1 + |a;|)2ll{|,|<^j(^^^-(:r) - r{x))h^ + 2 ^'^^f > 
<2de^ + 211(1 + |x|)2-'=/2l{|,|<^}(l + \x\)'^/'i^^'%x) - r (:r))|U. + 2 ^^^'^^^ +f ""''^ 



<2de' + Ci?(4+'^-^-)/2lj(l + Ixl)nA^^^^(x) - r{x))X{' + 2 ^^^(^^+f 



Jlk-2 

N,e\ 
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Proceeding as in Step 2 for the middle term and noting that for e S (0, 1), E[Mfc(/^ + jJ.^''^)] < 
C(l + Mfe(/)), we obtain, with d = 3, 

/)) < 6.^ + C^^il + M,(/))V2 + c'-+0il. 

Choosing e = and R = N'^l^^, we get a bound in A/'-(2fe-4)/5fe desired. 

Step 4: k = 7. It suffices to use the bounds we proved for any k <7. □ 
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